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Chapter  1 


Literature  review 

1.1  Quantum  dots,  their  physical  properties 
and  applications 

Quantum  physics  is  progressing  rapidly.  Quantization  and  spin  were  discov¬ 
ered  through  optical  studies,  as  were  other  fundamental  atomic  properties. 
With  the  advent  of  the  laser,  physicists  learned  how  to  manipulate  atomic 
wavefunctions  by  applying  coherent  optical  fields.  More  discoveries  followed. 
Now,  at  the  beginning  of  the  new  century,  optical  techniques  are  being  used 
to  explore  a  new  scientific  frontier:  the  atom-like  entities  known  as  quantum 
dots. 

A  quantum  dot  is  a  nanocrystal  made  of  semiconductor  materials  that  are 
small  enough  to  exhibit  quantum  mechanical  properties,  as  shonnd  in  (hgl), 
semiconductor  structures  in  which  the  electron  wavefunction  is  confined  in 
all  three  dimensions  by  the  potential  energy  barriers  that  form  the  quantum 
dots  boundaries  [4,  3].  A  quantum  dots  electronic  response,  like  that  of  a 
single  atom,  is  manifest  in  its  discrete  energy  spectrum,  which  appears  when 
electron-hole  pairs  are  excited.  Although  the  wavefunction  of  a  quantum  dot 
electron,  and  its  corresponding  hole,  extends  over  many  thousands  of  lat¬ 
tice  atoms,  the  pair-termed  an  exciton-behaves  in  a  quantized  and  coherent 
fashion. 

The  coherence  is  relatively  easy  to  detect  and  control  optically  -  for 
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Figure  1.1:  AFM  micrographs  of:  (a)  1  //m  x  1  fim  surface  imaging  of  In  As 
quantum  dots  on  GaAs/InP ,  (inset)  a  single  InAs  quantum  dot. 

two  reasons.  First,  the  superposition  of  the  ground  and  excited  states  de¬ 
phases  more  slowly  in  quantum  dots  than  in  higher-dimensional  semiconduc¬ 
tor  structures.  Second,  quantum  dots  have  large  dipole  moments  (50-100 
times  larger  than  those  of  atoms).  Thanks  to  these  advantages,  it  is  possible 
to  probe  and  manipulate  the  wavefunction  of  a  single  quantum  [9]. 

How  can  one  understand  the  excited  electronic  states  of  a  nanometer  sized 
semiconductor  crystallite,  given  that  the  crystallite  structure  is  simply  that 
of  an  excised  fragment  of  the  bulk  lattice?  This  question  is  motivated  by 
recent  experiments  on  chemically  synthesized  ’’quantum  crystallites,”  some¬ 
times  called  ’’quantum  dots,”  in  which  it  is  observed  that  the  optical  spectra 
are  quite  sensitive  to  size.  For  example,  bulk  crystalline  CdSe  is  a  semicon¬ 
ductor  with  an  optical  bandgap  at  690  nm ,  and  continuous  optical  absorption 
at  shorter  wavelengths.  However,  3540— diameter  CdSe  crystallites  contain¬ 
ing  some  1500  atoms  exhibit  a  series  of  discrete  excited  states  with  a  lowest 
excited  state  at  530nm  (1-3).  With  increasing  size,  these  states  shift  red  and 
merge  to  form  the  optical  absorption  of  the  bulk  crystal.  Electron  microscopy 
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and  BraggX-ray  scattering  measurements  show  that  these  crystallites  have 
the  same  structure  and  unit  cell  as  the  bulk  semiconductor.  Such  changes 
have  now  been  observed  in  the  spectra  of  many  different  semiconductors. This 
phenomenon  is  a  ”  quantum  size  effect”  related  to  the  development  of  the 
band  structure  with  increasing  crystallite  size.  Smaller  crystallites  behave 
like  large  molecules  (e.g.  polycyclic  aromatic  hydrocarbons)  their  spectro¬ 
scopic  and  photophysical  properties.  They  are  true  ”  clusters”  that  do  not 
exhibit  bulk  semiconductor  electronic  properties. 

Quantum  dots  possess  another  attractive  property.  Their  size,  shape, 
and  composition  can  all  be  ailored  to  create  a  variety  of  desired  properties. 
These  ”  artificial  atoms”  can,  in  turn,  be  positioned  and  assembled  into  com¬ 
plexes  that  serve  as  new  materials.  Many  types  of  quantum  dot  materials 
have  already  been  developed.  In  the  early  1980s  that  papers  in  the  Russian 
literature,  by  Ekimov  and  Efros,  gave  the  first  experimental  and  theoreti¬ 
cal  description  of  3D  quantum  confinement  with  semiconductor  nanocrystals 
[10,  11]. 

There  are  discovered  quantum  dots  in  a  glass  matrix  (they  are  developed 
a  growth  technique  of  the  semiconductor  microcrystals  in  a  glassy  dielec¬ 
tric  matrix.  This  technique  permits  to  vary  the  size  of  the  grown  micro¬ 
crystals  in  a  controlled  manner  from  some  tens  to  thousands  of  angstroms) 
and  in  colloidal  solutions  by  Louis  E.Brus  [12],  The  term  ’’quantum  dot” 
was  coined  by  M.  Reed.  [6]  (The  approach  he  used  to  produce  quantum 
dot  nanostructures  suitable  for  electronic  transport  studies  was  to  confine 
resonant-tunneling  heterostructures  laterally  with  a  fabrication-imposed  po¬ 
tential.  This  approach  embeds  a  quasibound  quantum  dot  between  two  quan¬ 
tum  wire  contact.)  Among  them  are  quantum  dots  formed  by  electrostatic 
gates.  Researchers  study  these  quantum  dots  with  far-infrared  spectroscopy 
and  with  transport  techniques,  such  as  conductance  measurements.  Great 
advances  have  been  made  with  electrostatic  quantum  dots,  but  this  type  of 
quantum  dot  does  not  interact  strongly  with  light  because  the  electrostatic 
gates  separate  the  electrons  and  holes,  severely  reducing  the  quantum  dots 
dipole  moment. 

Fortunately,  some  types  of  quantum  dots  do  interact  strongly  with  light 
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Core  —  CdSe,  CdS 
Shell  —  ZnS,CdS,ZnSe 
Amphiphilic  surface 

Cd  #Se/S 
Zn/Cd  #S/Se 

Figure  1.2:  Colloidal  quantum  dots  general  schematic 

and  are  being  studied  in  detail,  even  at  the  single  quantum  dot  level.  Grown 
as  part  of  a  larger  crystalline  structure, epitaxial  quantum  dots  form  spon¬ 
taneously  in  molecular  beam  epitaxy  [13].  Epitaxial  quantum  dots  are  es¬ 
pecially  attractive  to  researchers  because  they  lack  the  defects  caused  by 
etching,  regrowth,  and  other  processing  steps.  Another  advantage  of  epitax¬ 
ial  quantum  dots  is  that  they  lie  far  beneath  the  surface  of  the  surround¬ 
ing  material  whose  associated  states  are  potentially  meddlesome.  Colloidal 
chemistry  provides  yet  another  way  of  growing  quantum  dots. 

Colloidal  quantum  dots  -  semiconductor  nanocrystals  with  a  size  in  the 
range  10  -  20  nm  consisting  of  103-105  atoms  that  are  based  on  inorganic 
semiconductor  materials,  covered  with  a  monolayer  of  the  stabilizer  (’’coat”  of 
organic  molecules).  Colloidal  quantum  dots  combined  physical  and  chemical 
properties  of  molecules  with  optoelectronic  properties  of  semiconductors  (hg2 

)■ 

Colloidal  semiconductor  nanocrystals  are  synthesized  from  precursor  com¬ 
pounds  dissolved  in  solutions,  much  like  traditional  chemical  processes.  The 
synthesis  of  colloidal  quantum  dots  is  done  by  using  precursors,  organic 
surfactants,  and  solvents.  Heating  the  solution  at  high  temperature,  the 
precursors  decompose  forming  monomers  which  then  nucleate  and  generate 
nanocrystals.  The  temperature  during  the  synthetic  process  is  a  critical  fac¬ 
tor  in  determining  optimal  conditions  for  the  nanocrystal  growth.  It  must 
be  high  enough  to  allow  for  rearrangement  and  annealing  of  atoms  during 
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the  synthesis  process  while  being  low  enough  to  promote  crystal  growth.  The 
concentration  of  monomers  is  another  critical  factor  that  has  to  be  stringently 
controlled  during  nanocrystal  growth.  The  growth  process  of  nanocrystals 
can  occur  in  two  different  regimes,  ’’focusing”  and  ” defocusing” .  At  high 
monomer  concentrations,  the  critical  size  (the  size  where  nanocrystals  nei¬ 
ther  grow  nor  shrink)  is  relatively  small,  resulting  in  growth  of  nearly  all 
particles.  In  this  regime,  smaller  particles  grow  faster  than  large  ones  (since 
larger  crystals  need  more  atoms  to  grow  than  small  crystals)  resulting  in 
’’focusing”  of  the  size  distribution  to  yield  nearly  mono-disperse  particles. 
The  size  focusing  is  optimal  when  the  monomer  concentration  is  kept  such 
that  the  average  nanocrystal  size  present  is  always  slightly  larger  than  the 
critical  size.  Over  time,  the  monomer  concentration  diminishes,  the  criti¬ 
cal  size  becomes  larger  than  the  average  size  present,  and  the  distribution 
”  defocuses” . 

There  are  colloidal  methods  to  produce  many  different  semiconductors. 
Typical  dots  are  made  of  binary  alloys  such  as  cadmium  selcnide,  cadmium 
sulfide,  indium  arsenide,  and  indium  phosphide.  Dots  may  also  be  made 
from  ternary  alloys  such  as  cadmium  selcnide  sulfide.  These  quantum  dots 
can  contain  as  few  as  100  to  100,000  atoms  within  the  quantum  dot  volume, 
with  a  diameter  of  10  to  50  atoms.  This  corresponds  to  about  2-10  nm,  and 
at  10  nm  in  diameter,  nearly  3  million  quantum  dots  could  be  lined  up  end 
to  end  and  fit  within  the  width  of  a  human  thumb. 

Large  batches  of  quantum  dots  may  be  synthesized  via  colloidal  synthesis. 
Due  to  this  scalability  and  the  convenience  of  benchtop  conditions,  colloidal 
synthetic  methods  are  promising  for  commercial  applications.  It  is  acknowl¬ 
edged  to  be  the  least  toxic  of  all  the  different  forms  of  synthesis. 

In  a  semiconductor  crystallite  whose  diameter  is  smaller  than  the  size 
of  its  exciton  Bohr  radius,  the  excitons  are  squeezed,  leading  to  quantum 
confinement.  The  energy  levels  can  then  be  modeled  using  the  particle  in  a 
box  model  in  which  the  energy  of  different  states  is  dependent  on  the  length 
of  the  box.  Quantum  dots  are  said  to  be  in  the  weak  confinement  regime 
if  their  radii  are  on  the  order  of  the  exciton  Bohr  radius;  quantum  dots  are 
said  to  be  in  the  ’strong  confinement  regime’  if  their  radii  are  smaller  than 
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Exciton 

(electron-hole  pair) 


(b)  Band  gap 


Zero  point  vibrational 
(c)  energy  of  the  excited 
electron 


Zero  point  vibrational 
energy  of  the  hole 


Figure  1.3:  Exciton  energy  levels 

the  exciton  Bohr  radius.  If  the  size  of  the  quantum  dot  is  small  enough 
that  the  quantum  confinement  effects  dominate  (typically  less  than  10mn), 
the  electronic  and  optical  properties  are  highly  tunable  (fig4)  [reference]. 
Fluorescence  occurs  when  an  excited  electron  relaxes  to  the  ground  state 
and  combines  with  the  hole.  In  a  simplified  model,  the  energy  of  the  emitted 
photon  can  be  understood  as  the  sum  of  the  band  gap  energy  between  the 
occupied  level  and  the  unoccupied  energy  level,  the  confinement  energies  of 
the  hole  and  the  excited  electron,  and  the  bound  energy  of  the  exciton  (the 
electron-hole  pair). 

The  band  gap  can  become  larger  in  the  strong  conEnement  regime  where 
the  size  of  the  quantum  dot  is  smaller  than  the  Exciton  Bohr  radius  a*,*  as 
the  energy  levels  split  up.  where  a&* is  the  Bohr  radius  Tb  =  0.053  nm,  m  is 
the  mass. This  results  in  the  increase  in  the  total  emission  energy  (the  sum  of 
the  energy  levels  in  the  smaller  band  gaps  in  the  strong  confinement  regime 
is  larger  than  the  energy  levels  in  the  band  gaps  of  the  original  levels  in  the 
weak  confinement  regime)  and  the  emission  at  various  wavelengths;  which  is 
precisely  what  happens  in  the  sun,  where  the  quantum  confinement  effects 
are  completely  dominant  and  the  energy  levels  split  up  to  the  degree  that 
the  energy  spectrum  is  almost  continuous,  thus  emitting  white  light. 

The  exciton  entity  can  be  modeled  using  the  particle  in  the  box.  The 
electron  and  the  hole  can  be  seen  as  hydrogen  in  the  Bohr  model  with  the 
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Figure  1.4:  Splitting  of  energy  levels  for  small  quantum  dots  due  to  the 
quantum  confinement  effect.  The  horizontal  axis  is  the  radius,  or  the  size,  of 
the  quantum  dots  and  a*,*  is  the  Exciton  Bohr  radius. 


hydrogen  nucleus  replaced  by  the  hole  of  positive  charge  and  negative  elec¬ 
tron  mass.  Then  the  energy  levels  of  the  exciton  can  be  represented  as  the 
solution  to  the  particle  in  a  box  at  the  ground  level  (n  =  1)  with  the  mass 
replaced  by  the  reduced  mass.  Thus  by  varying  the  size  of  the  quantum  dot, 
the  confinement  energy  of  the  exciton  can  be  controlled.  There  is  Coulomb 
attraction  between  the  negatively  charged  electron  and  the  positively  charged 
hole.  The  negative  energy  involved  in  the  attraction  is  proportional  to  Ry¬ 
dberg’s  energy  and  inversely  proportional  to  square  of  the  size-dependent 
dielectric  constant  [15]  of  the  semiconductor.  When  the  size  of  the  semi¬ 
conductor  crystal  is  smaller  than  the  Exciton  Bohr  radius,  the  Coulomb 
interaction  must  be  modified  to  fit  the  situation. 

Although  the  above  equations  were  derived  using  simplifying  assump¬ 
tions,  the  implications  are  clear;  the  energy  of  the  quantum  dots  is  dependent 
on  their  size  due  to  the  quantum  confinement  effects,  which  dominate  below 
the  critical  size  leading  to  changes  in  the  optical  properties.  This  effect  of 
quantum  confinement  on  the  quantum  dots  has  been  experimentally  veriEed 
[16]  and  is  a  key  feature  of  many  emerging  electronic  structures  [17].  Besides 
confinement  in  all  three  dimensions  (i.e.,  a  quantum  dot),  other  quantum  con- 
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Figure  1.5:  Colloidal  quantum  dots  general  schematics. 


fined  semiconductors  include:  Quantum  wires,  which  confine  electrons  or  holes 
in  two  spatial  dimensions  and  allow  free  propagation  in  the  third. Quantum 
wells,  which  confine  electrons  or  holes  in  one  dimension  and  allow  free  prop¬ 
agation  in  two  dimensions. 

Quantum  size  effects  play  a  key  role  in  the  optoelectronic  properties  of 
quantum  dots.  The  energy  spectrum  of  the  quantum  dot  is  fundamentally 
different  from  the  bulk  semiconductor  .  Electron  in  a  nanocrystal  behaves 
like  a  three  -dimensional  potential  well.  There  are  several  stationary  energy 
levels  for  electron  and  a  hole  with  a  characteristic  distance  between  them. 
The  energy  spectrum  of  the  quantum  dot  depends  on  its  size  .  Similarly, 
the  transition  between  the  transition  energy  in  the  atom  ,  the  transition  of 
charge  carriers  between  the  energy  levels  in  the  quantum  dot  can  be  emitted 
or  absorbed  photon.  Transition  frequencies  ,  i.e  long  wave  absorption  or 
luminescence  easily  manage  changing  the  dimensions  of  the  quantum  dots. 
Therefore,  quantum  dots  are  sometimes  called  ’’artificial  points  .”  In  terms  of 
semiconductor  materials  that  can  be  called  the  ability  to  control  the  effective 
band  gap. 

There  is  another  fundamental  property  distinguishing  colloidal  quantum 
dots  of  semiconductor  materials  -  the  possible  existence  in  the  form  of  solu¬ 
tions.  (Fig. 5).  Synthesis  of  colloidal  quantum  dots  is  wide  as  possible  to  ob¬ 
tain  the  quantum  dots  on  the  basis  of  various  semiconductors  and  the  prepa¬ 
ration  of  quantum  dots  with  different  geometry  (shape).  Colloidal  quantum 
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dots  are  characterized  by:  the  composition,  size  and  shape.  Exsist  many 
forms  of  quantum  dots  spherical,  ellipsoidal,  nanocrystals  with  complex  ge¬ 
ometry,  multi-quantum  dots  [9]. 

The  various  types  of  quantum  dots  each  have  their  advantages  and  dis¬ 
advantages  and  allhave  served  as  useful  models  in  research  aimed  at  gaining 
a  better  understanding  of  theoretical  properties  of  quantum  dots. 

Researchers  have  studied  applications  for  quantum  dots  in  transistors, 
solar  cells,  LEDs,  and  diode  lasers.  They  have  also  investigated  quantum  dots 
as  agents  for  medical  imaging  and  as  possible  qubits  in  quantum  computing. 
The  Erst  commercial  release  of  a  product  utilizing  quantum  dots  was  the 
Sony  XBR  X900A  series  of  flat  panel  televisions  released  in  2013. 

Quantum  dots  are  particularly  significant  for  optical  applications  due  to 
their  high  extinction  coefficient.  In  electronic  applications  they  have  been 
proven  to  operate  like  a  single  electron  transistor  and  show  the  Coulomb 
blockade  effect.  Quantum  dots  have  also  been  suggested  as  implementations 
of  qubits  for  quantum  information  processing. 

The  ability  to  tune  the  size  of  quantum  dots  is  advantageous  for  many 
applications.  For  instance,  larger  quantum  dots  have  a  greater  spectrum-shift 
towards  red  compared  to  smaller  dots,  and  exhibit  less  pronounced  quantum 
properties.  Conversely,  the  smaller  particles  allow  one  to  take  advantage  of 
more  subtle  quantum  effects. 

Being  zero-dimensional,  quantum  dots  have  a  sharper  density  of  states 
than  higher-dimensional  structures.  As  a  result,  they  have  superior  trans¬ 
port  and  optical  properties,  and  are  being  researched  for  use  in  diode  lasers, 
amplifiers,  and  biological  sensors.  Quantum  dots  may  be  excited  within  a  lo¬ 
cally  enhanced  electromagnetic  held  produced  by  gold  nanoparticles,  which 
can  then  be  observed  from  the  surface  plasmon  resonance  in  the  photolu- 
minescent  excitation  spectrum  of  ( CdSe )  ZnS  nanocrystals.  High-quality 
quantum  dots  are  well  suited  for  optical  encoding  and  multiplexing  applica¬ 
tions  due  to  their  broad  excitation  profiles  and  narrow/symmetric  emission 
spectra.  The  new  generations  of  quantum  dots  have  far-reaching  potential 
for  the  study  of  intracellular  processes  at  the  single-molecule  level,  high- 
resolution  cellular  imaging,  long-term  in  vivo  observation  of  cell  trafficking, 
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tumor  targeting,  and  diagnostics. 

Quantum  dot  technology  is  one  of  the  most  promising  candidates  for 
use  in  solid-state  quantum  computation.  By  applying  small  voltages  to  the 
leads,  the  flow  of  electrons  through  the  quantum  dot  can  be  controlled  and 
thereby  precise  measurements  of  the  spin  and  other  properties  therein  can 
be  made.  With  several  entangled  quantum  dots,  or  qubits,  plus  a  way  of 
performing  operations,  quantum  calculations  and  the  computers  that  would 
perform  them  might  be  possible. 

Quantum  dots  may  be  able  to  increase  the  efficiency  and  reduce  the  cost 
of  today’s  typical  silicon  photovoltaic  cells.  According  to  an  experimental 
proof  from  2004,  quantum  dots  of  lead  selenide  can  produce  more  than  one 
exciton  from  one  high  energy  photon  via  the  process  of  carrier  multiplication 
or  multiple  exciton  generation  (MEG).  This  compares  favorably  to  today’s 
photovoltaic  cells  which  can  only  manage  one  exciton  per  high-energy  photon, 
with  high  kinetic  energy  carriers  losing  their  energy  as  heat.  Quantum  dot 
photovoltaics  would  theoretically  be  cheaper  to  manufacture,  as  they  can  be 
made  ’’using  simple  chemical  reactions.” 

There  are  several  inquiries  into  using  quantum  dots  as  light-emitting 
diodes  to  make  displays  and  other  light  sources,  such  as  ’’QD-LED”  dis¬ 
plays,  and  ’’QD-WLED”  (White  LED).  In  June  2006,  QD  Vision  announced 
technical  success  in  making  a  proof-of-concept  quantum  dot  display  and  show 
a  bright  emission  in  the  visible  and  near  infra-red  region  of  the  spectrum. 
Quantum  dots  are  valued  for  displays,  because  they  emit  light  in  very  specific 
gaussian  distributions.  This  can  result  in  a  display  that  more  accurately  ren¬ 
ders  the  colors  that  the  human  eye  can  perceive.  Quantum  dots  also  require 
very  little  power  since  they  are  not  color  filtered.  Additionally,  since  the  dis¬ 
covery  of  ’’white- light  emitting”  QD,  general  solid-state  lighting  applications 
appear  closer  than  ever.  A  color  liquid  crystal  display  (LCD),  for  example, 
is  usually  backlit  by  fluorescent  lamps  (CCFLs)  or  conventional  white  LEDs 
that  are  color  filtered  to  produce  red,  green,  and  blue  pixels.  A  better  solu¬ 
tion  is  using  a  conventional  blue-emitting  LED  as  light  source  and  converting 
part  of  the  emitted  light  into  pure  green  and  red  light  by  the  appropriate 
quantum  dots  placed  in  front  of  the  blue  LED.  This  type  of  white  light  as 
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backlight  of  an  LCD  panel  allows  for  the  best  color  gamut  at  lower  cost  than 
a  RGB  LED  combination  using  three  LEDs. 

Quantum  dot  displays  that  intrinsically  produce  monochromatic  light  can 
be  more  efficient,  since  more  of  the  light  produced  reaches  the  eye.QD-LEDs 
can  be  fabricated  on  a  silicon  substrate,  which  allows  integration  of  light 
sources  onto  silicon-based  integrated  circuits  or  microelectromechanical  sys¬ 
tems.  A  QD-LED  integrated  at  a  scanning  microscopy  tip  was  used  to  demon¬ 
strate  fluorescence  near-held  scanning  optical  microscopy  (NSOM)  imaging. 

Quantum  dot  photodetectors  (QDPs)  can  be  fabricated  either  via  solution¬ 
processing,  or  from  conventional  single-crystalline  semiconductors.  Conven¬ 
tional  single-crystalline  semiconductor  QDPs  are  precluded  from  integration 
with  flexible  organic  electronics  due  to  the  incompatibility  of  their  growth 
conditions  with  the  process  windows  required  by  organic  semiconductors.  On 
the  other  hand,  solution-processed  QDPs  can  be  readily  integrated  with  an 
almost  infinite  variety  of  substrates,  and  also  postprocessed  atop  other  inte¬ 
grated  circuits.  Such  colloidal  QDPs  have  potential  applications  in  surveil¬ 
lance,  machine  vision,  industrial  inspection,  spectroscopy,  and  fluorescent 
biomedical  imaging. 

A  variety  of  theoretical  frameworks  exist  to  model  optical,  electronic,  and 
structural  properties  of  quantum  dots.  These  may  be  broadly  divided  into 
quantum  mechanical,  semiclassical,  and  classical. 

1.2  Quantum  bound  states  in  classically  un¬ 
bound  regions 

Advances  in  fine-line  lithography  allow  creation  of  devices  which  conduct 
along  two-dimensional  surfaces,  or  quantum  wires,  with  nanoscale  width, 
which  allow  electrons  to  propagate  in  the  channels  formed  by  these  surfaces, 
but  require  the  electron  wave  function  to  vanish  on  the  boundary  of  the 
surface  [1,  27].  Quantum  wires  have  been  used  extensively  to  investigate 
quantum  interference  effects  [28,  29,  30,  31,  32,  33].  Quantum- mechanical 
coherence  in  disordered  metals  and  the  associated  interference  phenomena  are 
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at  the  heart  of  many  topics  which  have  generated  intense  interest,  including 
localization,  resonant  tunneling,  and  universal  conductance  fluctuations  [29] . 
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Figure  1.6:  Schematic  diagram  of  the  potential  formed  by  the  crossed  wires. 

Schult  et.  al  reported  in  Ref.  [1]  calculations  designed  to  determine 
whether  an  electron  can  be  trapped  at  the  intersection  of  the  perpendicularly 
aligned  quantum  wires,  due  to  a  quantum  effects.  The  presence  of  such 
trapped  electrons  could  modify  considerably  the  operation  of  devices  made 
of  these  wires.  It  should  be  noted  that,  as  such  surfaces  have  no  “classically 
forbiden”  region  (a  classical  particle  could  roll  freely  through  such  a  system), 
the  discovery  by  Schult  et.  al  [1]  (and  earlier  by  Lenz  et.  al  [34])  that 
such  systems  possess  a  bound  state  was  rather  surprising.  It  shows  the 
quantum  binding  of  a  classically  unbound  system,  and  thus  complements  the 
phenomenon  of  the  decay  of  a  classically  bound  stte  by  quantum  tunneling. 
Goldstone  and  Jaffe  [35,  36]  (and  Exner  [37,  38])  then  proved  the  remarkable 
result  that  at  least  one  bound  state  exists  for  all  two-dimensional  surfaces  of 
constant  width  (except  surfaces  of  constant  curvature,  which  have  no  bound 
state). 

The  authors  of  [1]  asserted  that  the  detailed  shape  of  the  potential  which 
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confines  the  electron  to  the  channel  is  not  important  in  determining  the 
qualitative  fact  of  the  presence  of  a  bound  state  at  the  intersection.  For  an 
initial  wave  function  which  is  odd  about  the  x  and  y  center  lines  and  even 
about  the  x  =  y  line,  the  energy  is  E2  =  3.72 Et.  No  bound  states  are  found 
which  are  od  about  the  x  =  y  line.  The  main  disadvantage  of  this  mesh- 
point  method  is,  that  it  gives  only  the  lowest  state  of  each  symmetry  class. 
Possible  higher-energy  bound  states  would  be  inaccessible.  For  these  the 
Authors  of  Ref.  [1]  reported  their  second  method  which  uses  an  expansion 
in  a  complete  set  of  solutions  of  the  differential  equation  AT  =  —  ArT  (for  a 
trial  value  of  k2)  in  each  of  the  five  large  rectangular  regions  shown  in  Fig. 
(2.1).  The  coefficients  in  this  expansion  should  be  chosen  to  match  T  and  its 
derivatives  at  the  boundaries  between  these  regions.  One  of  the  advantages 
of  this  method  is  that  it  can  be  extended  to  energies  above  the  propagation 
threshold.  Also  it  can  be  easily  extended  to  any  geometry  composed  of 
rectangles,  open  or  closed,  with  equal  of  different  width.  For  a  tight-binding 
approach  to  such  problems,  see  [40].  This  method  may  be  used  to  investigate 
the  importance  of  the  sharp  corners  in  forming  the  bound  state,  and  the 
amount  of  bend  required  to  form  a  bound  state  of  the  bent  wire. 

As  is  well  known  (see  [41]),  a  “bulge”  in  a  two-dimensional  surface  can 
be  mapped  into  one  dimension;  the  transverse  bulge  then  appears  as  an  ef¬ 
fective  local  attraction,  which  in  one  dimension  always  produces  a  bound 
state.  However,  it  was  surprising  to  find  that  a  bend  produces  an  effective 
attraction  similar  to  a  bulge.  The  existence  of  these  bound  states  can  be 
understood  qualitatively.  In  Refs.  [42,  43]  Carini  et.  al  examined  the  prop¬ 
erties  of  systems  containing  a  single  bend.  The  Authors  examined  in  Ref. 
[43]  bent  waveguides  which  support  several  bound  states.  They  investigated 
this  question  from  both  theoretical  and  experimental  perspectives,  by  cal¬ 
culating  energies  and  holds  for  such  states  in  bent  waveguides  (such  as  in 
Fig.  (1.8),  constructing  such  waveguides,  and  comparing  their  experimental 
properties  with  the  theoretical  predictions.  In  [41]  Carini  et.  al  examine  the 
case  of  a  system  with  two  right-angle  bends  as  shown  in  Fig.  1.7.  The  width 
of  the  straight  sections  is  W  and  the  height  of  the  bend  is  H.  First  they 
consider  the  case  where  the  straight  sections  of  the  wire  are  infinitely  long. 
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The  wave  function  for  the  electron  satisfies  the  equation 

(V2  +  k2)ifi(x,  y)  =  0,  if)\s  =  0,  (1.1) 

where  the  wave  number  k  is  related  to  the  energy  E  by  k2  =  2 m*E/h2. 
In  either  straight  section  the  requirement  that  the  wave  function  vanish  on 
the  boundary,  and  separability  of  the  Hamiltonian,  forces  the  y  dependence 
of  the  wave  function  to  be  of  the  form  sm(nny/W)  for  integer  n.  This 
transverse  quantization  condition  produces  an  energy  threshold;  the  lowest 
energy  allowed  for  free  propagation  is  Ethr  =  (hi:)2 /(2m*W2).  Both  the 
extra  space  in  the  bend(s)  of  the  wire,  and  the  bending  itself,  produce  an 
effective  attraction  which  supports  electron  bound  state (s)  in  the  region  of 
the  bend(s).  The  wire  of  Fig.  1.7  will  have  one  or  two  bound  states,  which 
appear  as  isolated  states  with  energy  below  the  threshold  energy  Etj ir.  The 
bound  state  wave  functions  will  be  largest  in  the  vicinity  of  the  bend,  and 
fall  off  exponentially  with  distance  from  the  bend  region. 


Figure  1.7:  (From  Ref.  [41])  Simplified  model  for  bent  quantum  wire.  In¬ 
finitely  long  wire  with  two  right-angle  bends,  the  width  arbitrarily  normalized 
to  1  and  the  aspect  ratio  (height /width  ratio)  R.  (a)  the  aspect  ratio  R  >  2 
is  the  “quantum  bend  discontinuity”  case  of  Ref.  [48].  For  the  purpose  of 
calculations  the  wire  is  divided  into  three  sections  labeled  I,  II,  and  III,  re¬ 
spectively.  (b)  The  aspect  ratio  1  <  R  <  2  is  he  “quantum  bend  continuity” 
case,  the  wire  is  divided  into  sections  labeled  /  and  II,  respectively. 
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Thus,  bent  quantum  wires  are  examples  of  quantum  systems  whose  bound 
states  do  not  arise  from  the  “traditional”  picture,  where  a  binding  potential 
creates  classically  allowed  and  forbidden  regions.  Here  the  boundary  condi¬ 
tions  (vanishing  of  the  wave  function  on  the  boundaries  of  the  wire)  give  rise 
to  transverse  quantization  conditions  which  produce  a  minimum  threshold 
energy  for  continuum  solutions.  Localized  bends  or  bulges  in  these  wires 
then  produce  effective  local  attractive  forces  which  give  rise  to  bound  states. 
As  the  effective  attraction  in  a  bent  wire  is  (to  lowest  order)  proportional  to 
the  square  of  the  curvature  of  the  wire,  [35,  36]  the  magnitude  of  the  binding 
energy  increases  as  the  curvature  increases. 

The  analogy  of  bent  two-dimensional  systems  and  rectangular  waveguides 
is  considered  in  Refs.  [42,  43].  Suppose,  one  has  some  two-dimensional  curve 
a  in  the  xy  plane,  which  possesses  a  scalar  field  0  satisfying  the  Helmholtz 
(or  Schrodinger)  equation  [V2  +  k 2]0  =  0,  with  i/j\s  =  0  on  the  boundary  S 
of  the  curve,  then  one  can  produce  a  rectangular  waveguide  by  translating 
the  curve  a  normally  in  the  z  direction.  E  and  B  fields  can  be  constructed 
from  -0  as  follows: 

E(x,  y)  —  ikz*l>(x,  y),  B(x,  y)  =  —  z  x  V-0  (1.2) 

The  E  and  B  fields  of  this  equation  will  satisfy  Maxwell’s  equations  and 
boundary  conditions  for  TE  modes  in  the  waveguide  (see  e.g.  [44]),  where 
the  electron  wave  number  k  is  related  to  the  frequency  /  =  oj/2ti  by 

k  ;  27rv^/ 
c 

Therefore  if  there  exist  bound  states  of  the  Schrodinger  equation  for  particles 
moving  in  the  two-dimensional  curve  cr  (i.e.,  solutions  of  the  wave  equation 
below  the  minimum  energy  for  free  propagation  of  waves  in  the  wire),  there 
will  be  analogous  confined  TE  modes  (solutions  of  E  and  B  fields  with  fre¬ 
quencies  below  the  cutoff  frequency  for  the  waveguide),  and  the  confined  E 
field  of  Eq.  (1.2)  will  be  described  by  the  same  scalar  function  0  which  con¬ 
stitutes  the  bound  state  wave  function  of  the  electron  in  the  quantum  wire, 


(1.3) 
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in  Eq.  (1.1) 

Having  shown  correspondence  between  wave  functions  for  the  electrons 
moving  freely  in  a  quantum  wire,  and  the  electric  field  for  TE  modes  in  a 
rectangular  waveguide  (see  also  [37,  38,  27,  28,  29,  45,  42]),  Carini  et.  al 
[41]  constructed  waveguides  of  corresponding  shape,  and  demonstrated  the 
location  of  the  confined  states  by  measuring  the  ratio  of  reflected  to  incident 
power  R(f)  as  a  function  of  frequency  /  for  microwaves.  By  understanding 
the  properties  of  the  waveguides,  one  can  gain  information  about  the  quan¬ 
tum  channels  or  wires.  They  have  shown  that  theoretical  calculations  give 
very  good  agreement  with  experiments  of  confined  electromagnetic  fields  in 
bent  waveguides.  They  showed  that  the  number  of  confined  electric  held 
modes,  and  their  location,  are  determined  by  the  geometry  of  the  waveguide. 


Figure  1.8:  Sharply  bent  two-dimensional  surface. 

Experimental  studies  have  been  carried  out  for  quantum  wires  with  two 
bends  by  Wu  and  collaborators  [46,  47].  Wang  and  collaborators  [48]  carried 
out  theoretical  calculations  of  the  conductance  for  electrons  in  this  geometry, 
and  compared  it  to  these  experimental  results. 

Obviously,  the  considered  cross-shaped  region  received  application  in  the 
held  of  connected  wave  guides  [21].  Trapped  quantum  modes  in  complex 
three-dimensional  systems  were  addressed  in  [22],  The  eigenstates  of  an  elec¬ 
tron  in  an  infinite  quantum  waveguide  (e.g.,  a  bent  strip  or  a  twisted  tube) 
are  often  trapped  or  localized  in  a  bounded  region  that  prohibits  the  elec¬ 
tron  transmission  through  the  waveguide  at  the  corresponding  energies.  We 
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revisit  this  statement  for  resonators  with  long  but  finite  branches  that  is 
called  finite  waveguides.  Although  the  Laplace  operator  in  bounded  domains 
has  no  continuous  spectrum  and  all  eigenfunctions  have  finite  L2  norm,  the 
trapping  of  an  eigenfunction  can  be  understood  as  its  exponential  decay 
inside  the  branches.  We  describe  a  general  variational  formalism  for  detect¬ 
ing  trapped  modes  in  such  resonators.  For  finite  waveguides  with  general 
cylindrical  branches,  we  obtain  a  sufficient  condition  which  determines  the 
minimal  length  of  branches  for  getting  a  trapped  eigenmode.  Varying  the 
branch  lengths  may  switch  certain  eigenmodes  from  non-trapped  to  trapped 
or,  equivalently,  the  waveguide  state  from  conducting  to  insulating.  These 
concepts  are  illustrated  for  several  typical  waveguides  (L-shape,  bent  strip, 
crossing  of  two  strips,  etc.).  It  was  also  concluded  that  the  well-established 
theory  of  trapping  in  infinite  waveguides  may  be  incomplete  and  require  fur¬ 
ther  development  for  applications  to  finite-size  microscopic  quantum  devices. 

1.2.1  Connected  curved  waveguides,  mesoscopic  sys¬ 
tems 

One  could  consider  a  pair  of  parallel  straight  quantum  waveguides  coupled 
laterally  through  a  window  of  a  width  1  in  the  common  boundary.  The 
authors  of  [21]  showed  that  such  a  system  has  at  least  one  bound  state. 
We  find  the  corresponding  eigenvalues  and  eigenfunctions  numerically  using 
the  modematching  method,  and  discuss  their  behavior  in  several  situations. 
They  also  discussed  the  scattering  problem  in  this  setup,  in  particular,  the 
turbulent  behavior  of  the  probability  flow  associated  with  resonances.  The 
level  and  phaseshift  spacing  statistics  shows  that  in  distinction  to  closed 
pseudointegrable  billiards,  the  present  system  is  essentially  nonchaotic. 

A  strong  motivation  to  study  such  bound  states  and  related  resonance 
effects  comes  from  recent  developments  in  semiconductor  physics,  because 
they  can  be  used  as  models  of  electron  motion  in  socalled  quantum  wires, 
i.e.,  tiny  strips  of  a  very  pure  semiconductor  material,  and  similar  structures. 
Let  us  briefly  recall  key  features  of  such  systems.  Characteristic  properties  of 
the  semiconductor  microstructures  under  consideration  are  small  size,  typi- 
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cally  from  tens  to  hundreds  of  nm,  high  purity,  which  means  that  the  electron 
mean  free  path  can  be  a  few  /jrn  or  even  larger,  and  crystallic  structure.  In 
addition,  boundaries  of  the  microstructures  consist  usually  of  an  interface 
between  two  different  semiconductor  materials;  the  electron  wavefuction  are 
known  to  be  suppressed  there.  Behavior  of  an  electron  in  such  a  mesoscopic 
system  structure  is,  of  course,  governed  by  the  manybody  Schrodinger  equa¬ 
tion  describing  its  interaction  with  the  lattice  atoms  including  the  bound¬ 
ary,  external  fields,  and  possible  impurities.  The  mentioned  properties  allow, 
however,  to  adopt  several  simplifying  assumptions.  As  we  have  said  the  mean 
free  path  is  typically  two  or  three  orders  of  magnitude  greater  than  the  size 
of  the  structure;  hence  the  electron  motion  can  be  assumed  in  a  reasonable 
approximation  as  ballistic,  i.e.,  undisturbed  by  impurity  scattering. 

1.2.2  Bound  states  in  twisting  tubes 

Here  we  should  also  add  a  paper  by  J.  Goldstone  and  R.  L.  Jaffe  [23].  It  is 
easy  to  see  that  the  uncertainty  principle  can  pro-  vide  binding  in  quantum 
systems.  Consider,  for  example,  an  infinite  cylindrical  tube  with  a  bulge 
in  the  middle.  A  particle  constrained  to  move  in  this  tube  can  have  lower 
transverse  momentum  and  therefore  lower  transverse  energy  within  the  bulge, 
where  it  is  less  severely  confined.  A  bound  state  results  because  the  quantum 
particle  would  have  to  squeeze  into  the  cylindrical  section  of  the  tube  in  order 
to  propagate  away  to  infinity.  Such  a  system  a  straight  tube  with  a  slowly 
changing  radius  can  be  mapped  into  a  one-dimensional  problem  in  which  the 
changing  radius  appears  as  a  changing  potential.  A  bulge  corresponds  to 
attraction,  and  in  one  dimension  an  attraction,  no  matter  how  weak,  suffices 
to  produce  binding. 

The  object  of  this  paper  is  to  demonstrate  that  bends  do  as  well  as  bulges. 
It  has  been  proved  that  an  infinite  tube  with  a  constant  normal  cross  section 
always  has  a  bound  state  if  it  bends.  Apparently,  a  bend  provides  a  region  in 
which  the  quantum  particle  can  relax,  lowering  its  momentum  in  comparison 
with  straight  regions,  though  this  is  not  as  obvious  as  it  is  for  bulges. 
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1.2.3  Kirchhoffs  Rule  for  Quantum  Wires 

One  particle  quantum  scattering  theory  on  an  arbitrary  finite  graph  with  n 
open  ends  and  where  we  define  the  Hamiltonian  to  be  (minus)  the  Laplace 
operator  with  general  boundary  conditions  at  the  vertices.  This  results  in  a 
scattering  theory  with  n  channels  has  been  observed  in  [24],  The  correspond¬ 
ing  on-shell  S-matrix  formed  by  the  reflection  and  transmission  amplitudes 
for  incoming  plane  waves  of  a  negative  energy  explicitly  was  given  in  terms  of 
the  boundary  conditions  and  the  lengths  of  the  internal  lines.  It  appears  to 
be  unitary,  which  may  be  viewed  as  the  quantum  version  of  Kirchhoffs  law. 
The  system  exhibits  covariance  and  symmetry  properties.  It  is  symmetric 
if  the  boundary  conditions  are  real.  Also  there  is  a  duality  transformation 
on  the  set  of  boundary  conditions  and  the  lengths  of  the  internal  lines  such 
that  the  low  energy  behavior  of  one  theory  gives  the  high  energy  behavior  of 
the  transformed  theory.  Finally  we  provide  a  composition  rule  by  which  the 
on-shell  S-matrix  of  a  graph  is  factorizable  in  terms  of  the  S-matrices  of  its 
subgraphs.  All  proofs  only  use  known  facts  from  the  theory  of  self-adjoint 
extensions,  standard  linear  algebra,  complex  function  theory  and  elementary 
arguments  from  the  theory  of  Hermitian  symplectic  forms. 

1.3  Novel  numerical  methods  development 

Recently  this  problem  has  given  a  boost  to  the  development  of  a  number  of 
novel  numerical  methods  [25] ,  where  the  autohers  surveyed  various  numerical 
methods  for  finding  solutions  of  quantum  confined  states,  especially  consid¬ 
ering  states  in  two-dimensional  (2D)  tubes,  or  2D  surfaces  that  are  confined 
in  the  transverse  direction  but  are  unconfined  in  the  longitudinal  direction. 
We  first  review  existence  proofs  for  bound  states  in  long  2-D  tubes.  They 
also  reviewed  various  methods  for  finding  such  states  and  we  discuss  the 
significance  of  these  eigenstates. 

When  inventing  problems  for  introductory  courses  in  quantum  mechanics, 
it  is  common  for  instructors  to  avoid  the  complications  of  space-dependent 
potentials  by  dealing  with  motion  in  a  region  with  constant  (i.e.,  zero)  po- 
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tential,  bounded  by  hard  walls.  By  doing  so,  one  can  focus  on  the  effects 
of  the  boundary  conditions  on  the  solutions.  Thus,  the  student  typically 
begins  with  a  one-dimensional  box,  after  which  he/she  can  study  a  box  in 
two  or  three  dimensions.  As  we  will  show  in  this  article,  solutions  on  two- 
dimensional  surfaces  exhibit  a  number  of  interesting  quantum  properties, 
symmetries,  and  degeneracies.  In  this  paper,  we  consider  a  number  of  simple 
two-dimensional  (2D)  hard-wall  geometries,  which  can  be  straightforwardly 
solved  numerically  by  undergraduates,  and  we  give  examples  of  the  types  of 
solutions  that  can  be  demonstrated  for  such  systems.  One  example  of  a  2D 
system  that  has  been  widely  investigated  is  a  2D  region  where  a  particle  is 
confined  to  a  cavity.  Such  systems  have  a  direct  analog  with  electromagnetic 
waves  in  a  2D  cavity.  For  electrons  moving  in  a  2D  cavity  in  a  quantum 
heterostructure,  such  a  state  is  commonly  referred  to  as  a  quantum  dot.l  In 
this  paper,  our  emphasis  will  be  on  particles  moving  in  long  tubes  or  chan¬ 
nels,  where  the  particle  is  confined  in  the  transverse  direction  but  can  move 
freely  in  the  longitudinal  direction.  Such  states  for  electrons  in  quantum 
heterostructures  are  referred  to  as  quantum  wires  or  electron  waveguides”. 
It  has  been  relatively  recently  understood  that  under  a  rather  general  set 
of  conditions,  confined  modes  will  exist  in  these  2D  tubes,  despite  the  fact 
that  confined  modes  would  not  exist  for  classical  particles  in  such  a  system. 2 
As  we  will  show,  any  localized  bulge  or  bend  in  an  otherwise  straight  chan¬ 
nel  or  quantum  wire  will  give  rise  to  a  stationary  state  which  has  a  lower 
energy  than  the  threshold  energy  for  free  motion  down  the  channel;  the  cor¬ 
responding  wave  function  is  thus  localized.  Two  popular  textbooks  have 
posed  problems  dealing  with  these  states,  but  the  approaches  used  in  these 
problems  involved  variational  methods  or  limiting  cases,  which  are  sufficient 
to  demonstrate  the  existence  of  a  confined  mode  but  which  can  give  a  rather 
poor  value  for  the  energy.  In  this  paper  we  discuss  various  numerical  meth¬ 
ods  by  which  one  can  obtain  accurate  answers,  and  in  certain  cases  extremely 
precise  values,  for  the  eigenstates  in  simple  geometries. 
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Statement  of  the  problem 

2.1  S-type  skew  symmetric  cross- wire  system: 
motivation  and  reasoning 

With  the  advance  of  fine-line  lithography,  devices  can  now  be  made  which 
conduct  along  two-dimensional  surfaces  shaped  into  channels  the  width  of 
which  can  be  as  narrow  as  75  nm.  Configured  in  four  terminal  junctions, 
such  quantum  wires  are  used  to  investigate  quantum  interference  effects. 

We  introduce  a  skew-symmetric  quantum  wire  intersection,  as  shown  in 
Fig.  [2.1],  such  that  one  of  the  two  wires  is  straight,  whereas  the  other  one 
is  curved.  This  system  is  more  likely  to  posses  a  quantum  bound  state, 
compared  to  the  other  types  of  geometry. 

Since  for  the  configurations  we  consider  there  is  no  trapping  classically, 
the  presence  of  a  localized  quantum  mechanical  state  is  not  an  obvious  phe¬ 
nomenon.  It  shows  the  quantum  binding  of  a  classically  unbound  system, 
and  thus  complements  the  phenomenon  of  the  decay  of  a  classically  bound 
state  by  quantum  tunneling.  The  problem  is  distinct  from  open  geometries 
explored  previously,  such  as  that  where  the  conducting  region  is  contained  by 
hyperbolas  x2y 2  =  const.  Because  of  the  ” pinched-off”  arms  of  that  shape, 
there  are  an  infinity  of  classical  periodic  orbits,  so  that  the  presence  there  of 
quantum  bound  states  is  not  surprising.  We  return  to  this  point  later.  We 
believe  that  the  detailed  shape  of  the  potential  which  confines  the  electron 
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Figure  2.1:  Schematic  diagram  of  the  potential  formed  by  the  crossed  wires. 


to  the  channel  is  not  important  in  determining  the  qualitative  fact  of  the 
presence  of  a  bound  state  at  the  intersection.  We  take  a  potential  which  is 
zero  inside  the  channels  and  infinite  outside,  so  that  the  Hamiltonian  is 


H  = 


(2.1) 


inside  the  well,  with  the  boundary  condition  that  the  electron’s  wave  function 
goes  to  zero  on  the  sides  of  the  channel.  We  have  used  two  methods  to 
calculate  the  bound-state  energy  and  wave  function  and  describe  them  in 
the  following  paragraphs.  One  of  these  methods  is  also  suitable  for  extension 
to  the  propagating  case,  and  we  show  results  from  such  calculations. 

We  also  are  looking  for  a  different  type  of  solution  in  the  ”arms“  or 
’’sleeves”  away  from  the  central  square.  Once  the  influence  of  the  other 
arm  is  negligible,  the  problem  becomes  identical  to  one-dimensional  infinite 
potential  well  with  well  known  sin(7rx/ (2a)).  Thus  we  seek  the  solution  with  a 
parabolic  profile  which  provides  a  better  description  of  the  realistic  quantum 
state  compared  to  the  ”‘trianlge”  shape. 
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There  could  be  two  basic  methods  to  calculate  numerically  the  bound- 
state  energy  and  wave  function  and  describe  them  in  the  following  para¬ 
graphs.  One  of  these  methods  is  also  suitable  for  extension  to  the  propagating 
case,  and  we  show  results  from  such  calculations.  We  conclude  with  some 
remarks  on  other  configurations,  and  on  the  two-electron  case.  Schematic 
diagram  of  the  potential  formed  by  the  crossed  wires  is  the  following:  The 
hatched  area,  together  with  symmetry  statements,  is  sufficient  to  display  in¬ 
formation  about  the  bound-state  wave  functions  and  is  the  area  shown  in 
Fig.  2.  The  regions  I  through  V  are  used  in  our  function  expansion  method. 
The  origin  of  the  x  and  y  axes  is  taken  to  be  located  at  the  lower  left  corner 
of  region  V.  The  system  obviously  possesses  four-fold,  rotational  symmetry. 

The  wave  function  goes  to  zero  on  the  sides  of  the  channel  (see  Fig.  2.1). 
Two  methods  are  used  to  calculate  the  bound-state  energy  and  wave  function. 
The  first  model  is  a  mesh  point  method  in  which  the  Authors  replace  the 
Schrodinger  equation 

Hif  =  Eif  (2.2) 

by  a  difference  equation  for  the  wave  function  evaluated  on  a  rectangular 
mesh  of  points  in  the  plane  of  wires  with  discrete  evolution  in  a  “pseudotime” 
variable,  t  [39].  This  equation  is  the  discretized  version  of  the  differential 
equation 

5^/5t  =  AT,  (2.3) 

where  A  is  the  Laplacian  operator.  Then  the  difference  equation  is  given, 
which  was  iterated  until  T  becomes  an  eigenfunction  of  the  discretized  A 
operator.  When  the  initial  wave  function  is  taken  to  be  symmetric  about 
the  center  lines  in  the  x  and  y  directions  and  about  the  diagonal  line  x  —  y, 
the  eigenvalue  found  below  the  propagation  threshold  with  the  corresponding 
energy  E  =  0.66 Et,  where  Et  =  h2 (n / w)2 / (2m*)  is  the  threshold  energy,  w 
is  the  width  of  the  channel. 
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2.2  Variational  principle  and  formalism 

Variational  principle,  or  Rayleigh- Ritz  principle,  is  an  approximation  method 
of  calculations  of  the  one-particle  eigenvalues  (mostly,  the  ground  state),  dif¬ 
ferent  from  the  perturbation  method  [26]. 

Let  |4/n)  be  the  eigenstate  of  PL.  We  introduce  some  normalized  nearby 
state  4') 

l*>  =  l*n>  +  |<Stt>  (2.4) 

so  that  (4'1'L)  =  1.  This  means  that  if  |4/n)  is  a  complete  orthonormal  set  of 
the  eigenstates  of  PL,  then 

m  =  I*n>  +  \SV)  =  N\V)  +  Y,  (2.5) 

n^m 

for  some  small  numbers  em.  Let  £  be  the  expectation  value  of  PL  in  |T).  The 
main  idea  of  the  variational  method  is  that  the  change  of  £  due  to  the  em 
vanishes  in  the  first  order: 

-  (Tn|L/|Tn)  =  (^n\H\8^)  +  (5^\H\^n)  +  (5T|H|5T)  (2.6) 
The  sate  is  normalized,  (4/|'L)  =  1.  Therefore, 

(5T|T)  +  (T|hT)  =  <f(T|T)  =  0.  (2.7) 

Ground  state  is  the  most  important  case.  Let  (To)  -  ground  state  and  it 
is  not  degenerate.  Then  PL\P>)  —  £,°|4/o). 


(T|77|T)  =  YEn\(^n)\2  >  (2.8) 

n  n 


The  ground  state  is  a  true  minimum.  Our  solution  is  based  on  this  fact. 
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2.3  Existence  of  a  quantum  bound  state  in  a 
rectangular  crossed-wire  system 

Now  we  present  a  solution  to  a  simpler  considered  problem,  originally  sug¬ 
gested  by  Krishna  Rajagopal  in  1974.  The  ground  state  energy  is  below  the 
lowest  possible  propagating  state  energy,  therefore  the  existence  of  a  quan¬ 
tum  bound  state  is  confirmed.  We  demonstrate  the  existing  solution  of  the 
cross-shaped  infinite  region  based  on  [2], 


2.3.1  Previously  used  trial  function  by  K. Rajagopal 

We  start  with  Schrbdinger  equation  for  the  considered  geometry: 


h2  kd2^(x,  y)  d2^(x,y)\ 
2 m  \  dx2  dy 2  ) 


E^(x,  y) 


(2.9) 


Following  the  variable  separation  procedure,  we  look  for  the  solution  in  the 
following  form: 

V(x,y)  =  X(x)Y(y) .  (2.10) 

This  results  in  the  following  expression: 


Y 


«92X  <92Y 

X- 


dx2  dy 2 
being  split  into  two  separate  equations: 


2mE 


XY. 


(2.11) 


d2X 
dx 2 
d2  Y 

dy* 


(2.12) 

(2.13) 


with  k2x  +  k2  =  The  general  solution  to  Eq.2.13  is  Y (y)  =  C\  cos  kyy  + 
C2  sin  kyy,  and,  in  order  to  satisfy  the  boundary  conditions  y(±a)  =  0,  we 
need  to  assume  ky  =  n^-  with  its  minimal  value  ky  —  Therefore,  the 
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propagating  particle  energy  must  satisfy 


S  > 


2m 


(2.14) 


We  conclude  that  every  state  with  its  energy  eigenvalues  being  less  then  the 
critical  energy. 


Ecr 


7 CK2 
8  ma2 


(2.15) 


represent  a  bound  state. 


Here  we  address  a  trial  function  for  our  problem,  first  suggested  by  Kr¬ 
ishna  Rajagopal  in  1974  [2]: 


(l-p)e- 
(l  ~  a )  e~“lyl 
(l  -  l)  e~a\x\ 
0 


|x|  <  a  &&  \y\  <  a 
|x|  <  a  &&  \y\  >  a 
|x|  >  a  &&  \y\  <  a 
otherwise 


(2.16) 


Obviously,  the  function  is  continuous  and  drops  to  zero  in  the  boundaries, 
and  therefore  satisfies  all  the  requirements  for  a  wave  function.  However,  the 
first  derivative  of  this  trial  function  is  broken  along  certain  lines,  referred  to 

as  roof  lines. 


2.3.2  Solution  based  on  variational  principle 

Here  we  briefly  review  the  solution  for  the  ground  state  energy  corresponding 
to  the  previously  considered  K. Rajagopal  trial  function. 

We  note  that  generally  the  solution  consists  of  three  parts  -  normalization 
of  the  wave  function,  bulk  integration  as  well  as  the  roof  lines  integration. 
While  the  first  two  steps  are  rather  straightforward,  the  roof-line  part  requires 
special  integration  techniques,  involving  generalized  function  properties.  We 
are  going  to  develop  general  integration  method,  leading  to  the  roof-line 
integrals  for  both  Rajagopal  and  S-type  trial  functions.  This  method  will  be 
elaborated  in  details  in  Appendix  A.  As  for  now,  we  will  just  summarize  the 
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results,  as  it  was  presented  in  [2], 

First,  we  not  that  the  whole  interior  of  the  cross-shaped  region  consists  of 
8  pairs  of  distinct  regions,  triangular  1/8  symmetric  cut  of  the  central  square 
and  a  half  of  each  sleeve  (this  parts  will  be  referred  to  as  ”1“  and  ”2“  ).  The 
normalization  condition  could  be  now  presented  as: 

8(I1  +  I2)  =  1  (2.17) 


(2.18) 


(2.19) 


(2.20) 

(2.21) 


(2.22) 


CHAPTER  2.  STATEMENT  OF  THE  PROBLEM 
where  J\  vanishes: 


*-W(‘-3>( £+£>(■ -D ■“> 
0 


In  contrast,  J2  term  is  as  follows: 


J2  =  ^ 


l~vNv,a)  {^  +  w-[=0] 


x=a  y= 0 


(2.24) 


( 1  -  -e~lx^a\  dxdy  =  -A2r)t 
\  a  /  6 


We  add  up  all  the  results,  and  without  the  roof  lines  the  ground  state  energy 
is: 

—  e~27  (2.25) 


(H)---A2^e 
3  m 


Now  we  summarize  the  results  for  the  roof  line  integration. 

The  roof  lines  in  this  case  are  classified  as  follows: 

•  right  arm  central  symmety  line  (x8)  K\  ; 

•  central  square  symmetry  lines  -  /C 2  ; 

•  central  square  boundary  lines  -  1C  1  ; 

After  the  applying  generalized  funciton  properties,  the  integration  results  are 
as  follows: 


9  h2 

—  e~2^A2 

(2.26) 

m'y 

2  h2 

—  e~2^A2 
m 

(2.27) 

-  (37  -  1)  A2 

(2.28) 

Finally,  after  taking  into  account  all  the  required  components  and  the 
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proper  normalization,  the 


(V)  =  3ft2  72  +  27  +  3 

ma 2  6  +  II7 


(2.29) 


Consequently,  we  are  looking  to  the  variational  parameter  7  value,  minimiz¬ 
ing  the  ground  state  energy: 


d(H) 

dy 


=  0. 


(2.30) 


Here  we  are  obviously  looking  for  a  minimum,  the  maximum  does  not  exist. 
This  minimum  condition  (Eq.2.30)  results  in 


(6  +  ll7)(27  +  2)  =  ll(72  +  7  +  3)  (2.31) 


and 


7  = 


-6  +  V36  +  21 


11 


7  =  0.940 


and  the  corresponding  energy  is 


E0  =  (H)  = 


3ft2  3(7  +  1) 


ma “ 


11 


6  ft2  ft2 

-A,  (7+  1)  =  1-058— 


11  my2 


ma * 


(2.32) 


(2.33) 


This  result  is  obviously  an  approximation.  The  true  value  of  the  ground 
state  energy  is  less  than  Eo.  However,  the  obtained  energy  values  is  below 
the  threshold  energy  obtained  in  the  previous  section: 


-  ^  ft2 

8  ma 2 


1.2337- 


ft2 


ma z 


(2.34) 


Consequently,  the  existence  of  a  quantum  bound  state  is  proved.  As 
the  next  step,  we  will  present  a  novel  trial  wave  function,  leading  to  a  lower 
ground  state  energy  value. 


Chapter  3 

S-type  skew  symmetric  crossed 
wires:  ground  state 


In  this  Chapter,  the  calculations  are  discussed  to  determine  whether  quantum 
effects  can  cause  the  trapping  of  an  electron  at  the  intersection  of  S-type  skew 
symmetric  cross-wire  system.  The  sketch  of  the  system  is  shown  in  Fig. 3.1. 
Such  electron  trapping  could  modify  considerably  the  operation  of  devices 
made  of  these  wires  by  changing  their  transport  and  optical  properties,  such 
as  electron  mobility  and  localization  parameters,  that  are  important  for  op¬ 
tical  sensing  in  particular. 

The  principle  problem  here  is  to  develop  a  new  trial  function  that  would 
lead  to  a  lower  value  of  the  ground  state  energy,  and  then  to  calculate  the 
corresponding  bound  state  energy  of  the  system.  Obviously,  the  wave  func¬ 
tion  must  be  smooth  (the  function  itself  and  its  first  derivative  must  both  be 
continuous),  otherwise  quantum  current  can  not  be  defined  correctly.  The 
point  is  that  one  or  a  few  spatial  derivatives  (the  gradient  components)  of 
the  wave  function  is  an  essential  and  necessary  part  of  the  current  definition 
in  Quantum  Mechanics.  Consequently,  the  current  cannot  be  calculated  if 
the  wave  function  is  not  smooth  due  to  the  non-existence  of  the  wave  func¬ 
tion  first  derivatives  in  some  directions.  However  we  consider  a  stationary 
states  calculations  with  no  current  here  and  so  in  our  trial  wave  function 
there  could  be  lines  of  discontinuity  of  the  first  derivative.  The  latter  will 
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be  referred  to  as  ’’roof  lines”  in  what  follows.  Roof  lines,  or  lines  with  the 
discontinuous  spatial  derivative  of  the  wave  function,  are  the  central  point 
here.  The  system  behavior  along  these  lines  is  clearly  unphysical  since  the 
first  derivative  of  the  wave  function  does  not  exist.  In  principle  this  leads 
to  singularities  (non-standard  behavior)  of  a  quantum  current,  but  this  does 
not  affect  the  bound  state  calculations  with  no  current  we  are  interested  in. 
Also  such  regions  (or  lines)  are  normally  associated  with  5-function  (infinite) 
type  of  potential,  which  has  nothing  to  do  with  the  potential  profile  of  our 
problem  on  the  wave  function  spatial  domain. 

We  do  try  to  minimize  the  number  and  the  impact  of  the  roof  lines  in  our 
variational  model.  To  this  end,  we  are  looking  for  a  better  type  of  solution  in 
the  arms  (or  sleeves)  away  from  the  central  square  (see  Fig. 3.1),  where  the 
influence  of  the  other  arm  is  negligible  and  so,  the  problem  becomes  identical 
to  one-dimensional  infinite  potential  well  with  the  well-known  ground  state 
wave  function  solution  of  the  form  y/l/asin[a:/(2a)]  (2a  is  the  sleeve  width). 
The  simplest  wave  function  profile  to  mimic  this  behavior  is  the  parabolic 
profile  which  we  chose  over  here  to  provide  a  better  description  to  the  realistic 
quantum  state. 

3.1  The  trial  wave  function 

First,  we  specify  all  regions,  where  the  wave  function  is  specified  by  a  separate 
formula  (see  Fig. 3.1): 

Ri  (central  curved  rhombus) 

The  region  of  the  highest  probability  density  is  as  follows; 

(  (|x|  -  a)2  +  (\y\  -  a f  >  a 2 
<  |a:|  <  a 
{  \y\<a 


with  the  constant  wave  function 
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Rl  R2,{a,b,c,d} 
R441-8}  J 


R43 

R4,5 

exp.  decay 

-< 

exp.  d 

-  » 

Rf4  " 

R4,6 

profile 

potential 


Figure  3.1:  General  schematics  and  geometry  of  the  problem  with  finite  incli¬ 
nation  angle  a.  The  regions  of  the  different  wave  function  spatial  dependence 
are  specified  by  different  colors. 


Ti  (x,y)  =  Afe  7  =  const  , 


(3.1) 


Since  the  region  is  considerably  larger  compared  to  the  case  of  Rajagopal 
wave  function,  our  present  normalization  constant  is  expected  to  exceed  the 
previously  established  normalization: 


MRj 


9y 


,27 


2a2 lly  +  6 


(3.2) 


f?2  (four  quarter-circles) 

The  remaining  part  of  the  central  square  could  be  specified  as 


(|x|  —  a)2  +  (\y\  —  a)2  <  a 2 
|x|  <  a 

M  <  « 


(3.3) 


CHAPTER  3.  S-  TYPE  SKEW  SYMMETRIC  CROSSED  WIRES:  GRO  UND  STATE33 


Figure  3.2:  Roof  line,  representing  the  beginning  of  the  exponential  decay. 
The  wave  function  is  demonstrated  for  x  of  the  range  from  —0.92  a  to  0.92  a, 
and  y  -  from  0  to  2.24  a.  The  wave  function  is  presented  as  a  density  plot  using 
Temperature  colormap,  in  which  the  red  color  corresponds  to  the  highest 
probability  density  and  the  blue  one  -  to  the  lowest.  This  convention  is 
adopted  for  all  the  density  plots  in  the  present  Thesis. 


with  another  wave  function 


A  =  \/(a-  |i|)2  +  (a-  |t/|)2 


(3.4) 

(3.5) 


so  the  wave  function  depends  only  on  the  distance  from  the  adjoining  ver¬ 
tex  of  the  central  square.  Obviously,  the  function  is  isotropic  (no  angular 
dependence)  in  the  translated  coordinates. 


i?3  (two  circular  segments) 

•  The  upper  circular  segment  (C.S.)  is  located  in  the  following  region: 


V  >  a 

y  <  —Tx  +  a(l  +  T) 

(x  —  a)2  +  (y  —  a)2  <  (2a)2  , 


(3.6) 
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•  the  lower  one: 

V  <  ~a 

y  >  -Tx  -  a(  1  +  T)  (3.7) 

(x  +  a)2  +  {y  +  a)2  <  (2a)2 

where  T  =  tan  a  and  C  —  1/T  —  cot  cc.  The  corresponding  wave  function  is 
given  as: 

l-(02)e-  (3.8) 

A'  =  y/(a  =f  x)2  +  (a  =f  y)2  ,  (3.9) 

where  the  =F  sign  correspond  to  the  upper  and  lower  regions  respectively,  so 

that  the  corresponding  vertices  of  the  central  square  are  {a,  a}  and  {—a,  —a}. 
The  dependence  is  very  similar  to  that  of  i?2  region,  however  the  wave  func¬ 
tion  is  equal  to  zero  for  both  A'  =  0  and  A'  =  2 a  to  ensure  zero  probability 
density  at  all  boundaries  (on  the  both  sides). 


ty3(x,y)=Af  1 


R4  (Eight  semi-infinite  rectangles  -  ’’sleeves”) 

At  this  point  we  need  to  account  for  the  difference  between  the  horizontal 
(i?4,{3,4,5,6})  and  inclined  vertical  (-R4, {1,2, 7,8})  sleeves. 

The  region  R4  consists  of  the  lower  inclined  sleeve: 

{y  <  —  Tx  —  a(l  +  T) 

y  <  cx  —  a(l  —  c)  (3.10) 

y  >  cx  —  a(l  —  c  +  2vA  +  c2) 

and  the  upper  one: 

{y  >  —Tx  +  a(l  +  T) 

y  >  cx  +  a(l  —  c)  (3-11) 

y  <  cx  +  a(l  —  c  +  2y/l  +  c2) 


We  write  the  wave  function  in  such  way  that  the  equations  would  rep¬ 
resent  similar  coordinate  dependence  in  the  slanted  {£,  77}  local  coordinate 
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❖ 

0 

y=5.oo  Rj 

O 

<> 

0 

O 

Figure  3.3:  Spherical  trial  function  with  strong  attenuation. 


systems  for  each  sleeve,  presented  at  Fig.l: 


#4, {1,2}  =>■  ^4,/  =  Af 

—  cx  +  y  —  a(l  —  c  +  \/l  +  c2) 


-)1 

e 

a J  _ 

-7(1 +v/A 


V  = 


and,  in  a  similar  manner: 


VTTc? 

| Tx  +  y  -  a(l  +  T)\ 


W+r2 


Ri,{7,8}  =>  ^4,,/  =  A/ 

—  cx  +  y  +  a(l  —  c  +  y/l  +  c2)  | 


V] 

e 

a J  _ 

•7(1+^/“) 


e  = 


V  = 


Vi  +  c2 
\Tx  +  y  +  a(l  +  T)\ 


Vl  +  Ti 

In  the  case  of  a  symmetric  quantum  dot  (a  =  0): 


Ra,v  =>  ^4,v  =  AT  (  1  —  (  - 


x 


*-i\y\/a 


(3.12) 

(3.13) 

(3.14) 


(3.15) 

(3.16) 

(3.17) 

(3.18) 
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According  to  our  model,  only  inclined  vertical  sleeves  -R4>{  1,2, 7, 8}  alter  with 
angle  a ,  whereas  the  horizontal  sleeves  #4, {3,4, 5, 6}  are  not  modified  (as  well 
as  the  central  square): 


Ri,h  =>  'I U,h  =  N 


0—71*1/° 


(3.19) 


One  can  easily  verify  that  for  a  =  0  the  circular  segments  i?3  disappear 


Figure  3.4:  Density  plot  of  the  trial  function  for  finite  values  of  the  inclination 
angle  a.  The  plots  are  presented  for  three  different  angles  a  =  30°,  50° 
and  70°  and  for  two  values  of  the  decay  constant  -  7  =  0.4  and  7  =  0.7, 
corresponding  to  the  variational  minimum  of  the  Hamiltonian  expectation 
value. 


and  the  inclined  vertical  sleeves  wave  functions  (3. 12), (3. 15)  tend  to  the 
expression  of  the  wave  function  in  the  vertical  sleeve  (3.18).  These  wave 
functions  are  different  from  the  one  obtained  by  Ragajopal  as  predicted. 
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Figure  3.5:  Roof  line,  representing  the  beginning  of  the  exponential  decay. 
The  wave  function  is  demonstrated  for  x  of  the  range  from  —0.92  a  to  0.92  a, 
and  y  -  from  0  to  2.24  a. 

3.2  Variational  approach  to  the  expectation 
values 

In  order  to  obtain  the  new  normalization  coefficient,  we  perform  the  integra¬ 
tion  of  the  probability  density  |\Eh(:c,  y)\2- 

R[.  The  wave  function  i (re,  y )  is  constant  in  the  central  curved  rhombus. 
First,  we  need  to  obtain  the  area  of  R\ 

Ai  =  (2a)2  -  4  =  (4  -  vr)a2  -  0.858a2  (3.20) 

and  the  corresponding  probability  density: 

Ti  (x,y)  =  A/’e-7  =>■  |Ti(a;,r/)|2  =  J\f'2e~ 27  (3.21) 

Xx  —  f  dx  f  dy  |4>i(x,  y)\2  =Af2A1e~2'y  =  (4  —  7r)A/"2a2e~27  . 


i?2-  Region  i?2  consists  of  four  equivalent  quarter-circles.  We  demon¬ 
strate  the  result  for  the  quarter-circle  with  the  center  at  {a,  a},  so  that  the 
contribution  from  R2  region  is  expected  to  be  4X2.  Correspondingly,  we  find 
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Figure  3.6:  3D  plots  of  the  new  trial  wave  function. 


the  wave  function  integral: 


tt2(z,j/)=.A/r^l-^l--J  )  e-7 

a  ?r/2 


r\  ■ 


(3.22) 


To  =  AT2  J  rdr  J  dcj)  ^1  —  ^1  —  -  j  ^  e  27  =  J\f2^na2e  27  , 


0  0 
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where  r  is  the  distance  from  {a,  a}-vertex  of  the  central  square  to  the  point 
of  observation  for  our  choice  of  the  sample  region. 

Region  R->  vanishes  for  a  =  0. 

R4.  We  consider  a  sample  horizontal  sleeve  and  keep  in  mind  that  the 
actual  contribution  of  this  region  comes  as  4X4 


(3.23) 

e— 27 x/a  _ 


The  total  integral  now  could  be  written  as: 

1 1  39n2p-2')'" 

X  =  Tx  +  4Z2  +  4X4  =  A f2  (4  -  7T )a2e“27  +  —  na2e~^  +  — -  =  1 

15  15q 

r2  _  157  e27  _  3.75  7  e27 

~~  4a2  8  +  (15  -  tt)7  “  8  +  11. 867  a2  ' 

Obviously,  our  normalization  coefficient  is  expected  to  be  considerably  smaller 
than  NrT  For  7  =  0.94,  leading  to  the  variational  minimum  of  the  Hamilto¬ 
nian  expectation  value  for  y),  the  normalization  coefficient  would  be 

A/"2  =  1.207/a2  in  comparison  to  ( J\fRj )2  =  1.697/a2. 


3.2.1  Bulk  Integration 

Now  we  are  at  position  to  evaluate  the  expectation  value  of  the  Hamiltonian 
everywhere  inside  the  cross.  As  a  matter  of  fact,  we  need  to  account  only  for 
two  different  integrals  and  £>4,  corresponding  to  R2  and  i?4  regions.  The 
actual  expression  is  as  follows: 

h2  /  d'2  d2  \  h2 

w*  =  -ss  +  ^ (3’24> 

Important.  First,  we  note  that  the  bulk  integral  B\  =  0,  since  the  wave 
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function  is  constant  in  region  Ri  (central  rhombus).  In  spite  of  the  fact  that 
the  components  of  the  gradient  (spatial  derivatives)  are  not  equivalent,  the 
general  convex- concave  symmetry  is  not  broken.  In  the  case  of  Rajagopal 
trial  function  the  situation  is  similar. 

In  each  case  we  need  to  evaluate  the  following  quantity: 


Bt  =  JdxdyU*,y)(a^1  +  d2nx'y) 

Ri 


dx2 


dy 2 


(3.25) 


since  the  wave  function  is  always  real.  The  latter  statement  also  leads  to  the 
fact  that  any  quantum-mechanical  current  for  the  suggested  trial  function  is 
equal  to  zero  everywhere  in  the  system. 


£>2  -  quarter-circles 

Since  the  wave  function  ^(r)  has  only  radial  dependence,  it  is  advisable  to 
take  its  Laplacian  as  well  as  the  spatial  integral  in  polar  coordinates: 


V2T2(r) 


d2R(r)  l<9T(r) 


b2  = 


7t/2 

rdr  [  d(f)  V2T2(r)  T2(r) 


o  o 


— -A r\e-^ 

6 


(3.26) 


B4  -  semi-infinite  rectangles  (sleeves) 

First,  we  evaluate  the  Laplacian  of  the  trial  function,  considering  one  model 
upper  vertical  sleeve  with  y  >  a 


V2R4(x,y)  =  AfV2 


g— 72// a 


2  )a2  —  ary2 


72//a 


(3.27) 

Here  we  first  encounter  the  situation,  where  both  x—  and  y—  second  deriva¬ 
tives  are  non-zero,  namely,  while  the  d 2 /dy2  is  the  result  of  the  exponential 
decay  (derivatives  of  all  orders  are  non-zero);  d2 / dx2  reflects  the  parabolic 
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transversal  profile  of  the  wave  function.  The  situation  was  not  encountered 
with  the  Rajagopal  trial  function. 

Consequently,  the  integral  is  in  the  form  of 

a  oo 

B4=  I  dx  J  dy^4(x,y)V2^4(x,y)  =  AN2  ^  e~27  (3.28) 

—a  a 

Also  it  is  interesting  to  mention  that  here  in  contrast  to  tyR:>(x,y),  the  pure 
term  (without  normalization)  B4/J\f 2  — y  oo  for  7  — >•  0,  which  corresponds 
to  the  infinite  region  of  finite  concavity  of  the  sleeves.  The  singularity  is 
compensated  by  the  similar  behavior  of  the  wave  function  normalization. 

Now  we  can  summarize  the  results  of  the  bulk  integration: 

w = -'L  + 144 = -4  ^ + =<3’29> 

=  —  N2— - —  (I672  —  57T7  —  40)  e-27 
loym 

So  far  we  just  note  that  for  7^  =  0.94  we  would  obtain  (IT)  =  +0.440  A f  ft2 /m, 
whereas  for  ^^(x,  y)  this  term  is  purely  negative  for  any  value  of  7. 


3.2.2  Roof  Line  Integration 

One  of  the  most  important  features  of  the  new  trial  wave  function  is  minimal 
number  and  complexity  of  the  roof  lines.  In  fact,  there  is  only  one  type  of 
such  line  Ci,  associated  with  the  beginning  of  the  exponential  decay  in  the 
sleeves. 

We  recall  the  analytic  expressions  of  the  function  4/2 (x,y)  in  the  close 
vicinity  of  C\ 


-7 


T2(a:,  y)  =  Af  (l  -  (1  -  A/a)2)  e 
A  =  \/{a-  \x\ )2  +  (a  -  \y\ )2  , 


(3.30) 
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so  we  can  expand  it  in  terms  of  the  small  x— displacement  from  namely 
Sx  =  (a  —  x)  ja 

^x<a{x,y)=N  e~7  ,  (3.31) 

LV  J  a~y  J 

Obviously  the  expansion  is  not  valid  at  y  =  ±a.  The  physical  reason  to  that 
is  that  both  components  of  the  wave  function  gradient  (spatial  derivative) 
are  not  continuous  However,  the  vertices  of  the  central  square  do  not  give  any 
contribution  to  any  of  the  integrals  since  the  wave  function  \Eq(±a,  ±a)  =  0. 
Now  the  generic  wave  function  next  to  C\  is  presented  as: 

F0(x,y)  =  Ar  |(l  -  0  )  dx>a+  -  0  )  +  0*<a} 

xexp  -7  (ex<a  +  ^0x>a)  • 

Anzats.  Since  the  only  source  of  non-smooth  C°-type  behavior  is  the 
beginning  of  the  exponential  decay,  and  we  are  collecting  only  5-function 
terms,  the  sought  second  derivative  could  be  written: 


F0(x  ->•  a,  y)  =  M  ^1  -  0  ^  e  7 


The  function  is  smooth  (no  roof  lines)  for  7  — y  0,  so  we  can  write 


<92T 

dx2 


a 


S(x  —  a) 


(3.32) 
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*  >-(!)') -¥ 


(3.33) 


The  roof  line  integral  could  be  expressed  in  the  final  form: 


a  a+e 


jC  =  4  •  (  —  - —  )  J\T2  dy  dx^(x  a,  y) 


d24>(x  — >  a,  y) 


—a  a— e 


(3.34) 


£ ^  hyti-  (3.35) 

m  /  V  Va/  /  a  15  m 


3.3  Hamiltonian  Expectation  Values 

3.3.1  Zero  inclination  angle 

Now  we  find  the  Hamiltonian  average  as  the  sum  of  the  bulk  and  the  roof 
integrals: 


3  =  4  +  B2)  =  “I (16o/2  _  57r7  _  40) e  27  (3-36) 


32t2 

£  = 

15  m 


so  that  the  expectation  value: 


(3.37) 


h2  e-27 

(H)  =  B  +  C  =  —M2—-  (I672  +  57T7  +  40)  . 
m  157 


(3.38) 


Taking  into  account  the  normalization  constant 


U2  =  =4 


157  e27  3.75  7  e27 


4a2  8  +  (15  —  7r)7  8  +  11. 867  a2 


(3.39) 


<«>  = 


h2  I672  +  57T7  +  40 
ma 2  6O7  —  47T7  +  32 


we  finally  obtain 


(3.40) 
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Applying  the  variational  principle,  we  find  the  values  of  the  exponential  de¬ 
cay  constant  ’’gamma”,  corresponding  to  the  minimum  of  the  ground  state 
energy.  This  results  in  7  =  0.840  and  S  =  0.898 h2/(ma2). 


Figure  3.7:  Energy  expectation  value  as  a  function  of  the  exponential  decay 
7,  which  is  chosen  as  the  variational  parameter. 


3.3.2  Finite  non-zero  inclination  angle. 

As  we  already  mentioned,  the  number  and  the  type  of  the  roof  lines  are  not 
altered  in  the  case  of  y-asymmetric  quantum  dot.  Correspondingly,  we  need 
to  supplement  the  normalization  with  term  2  xl3,  and  bulk  integrals  -  with 
£b,  from  region  II, >: 


^3  fa,  y)  =  A/"  1 


A'  =  \/(a  =F  x)2  +  (a  =F  y)2  , 


(3.42) 
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so  that  the  integrals 


a  2a  2 a  2 

13  =  AT2  j  d(f>  J  rdr  |T3(r)|2  =  ae-27  J  ^1  —  ^1  —  -  j  ^  rdr$A3) 


o  o 

l3=Af^aa2e~^ 

15 


and 


V2T3(r)  = 


d 2T3(r)  l<9T3(r)  ,2  fa 


dr 2  +  r  dr  ^ a2  \r  2'& ^ 


2  a 


rdr 


1  -  f  1  -  - 

a 


4d--2]=-S 


a 2  vr 


a  a 


B3  =  /  rdr  /  <70  T2(r)  V2T2(r)  =  —JV2-ae  2,1 


o  o 


and,  finally, 

Now  we  summarize  all  the  terms, 


h2  4  h2  -2  _2 

—  -  fi3  =  +-  AT  a;  e 
2m  3  m 


X3  =  JV2 (4  —  7T )a2e  27 
X2  =  JV2— 7ra2e”27 


60 

,16 


X3  =  AA— aa“e 
15 


2q~27 


1a=JV2 


15y 


2  -27 
-a  e  7 


(3.45) 

(3.46) 

(3.47) 

(3.48) 

(3.49) 

(3.50) 

(3.51) 


so  that 


X  =  J^XA  =  1  (3.52) 

A 

1  4a2 

—j-  =  Xi  +  X2  +  X3  +  X4  =  — —  (8  +  7(15  —  7T  +  8a))  e  2 7  (3.53) 
AM  157 


The  normalization  constant  JV  is,  therefore,  equal  to: 
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Figure  3.8:  Ground  state  energy  for  the  skew  symmetric  cross-wired  system. 
The  upper  panel  shows  how  the  variational  minimum  for  the  decay  coefficient 
7  depends  on  the  inclination  angle  a  and  the  lower  two  plots  represent  the 
angular  dependence  of  the  ground  state  energy  and  the  bound  state  energy 
respectively. 
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JV  =  e 7  (8  +  7(15  —  7 r  +  8a))  1//2 
The  bulk  and  roof  line  integrals,  therefore,  give: 


(3.54) 


<«>  =  (5,r  +  8  (5“  +  ®  +  27))  e'27  (3-55) 

which  result  in  the  following  final  equation  to  determine  the  ground  state 


energy: 


(H}  = 


K2  57T7  +  8(5  +  5cry  +  2y2) 


(3.56) 


4  ma?  8  +  q(15  —  n  +  8cd) 

Now  we  are  looking  for  the  angular  dependence  of  the  variational  minimum 
of  the  ground  state  energy.  Since  we  chose  the  exponential  decay  coefficient 
7  to  be  that  variational  parameter,  its  angular  dependence  is  now  obtained 
by  simple  differentiation: 

=  0  (3.57) 


This  results  in  the  following  angular  dependence  of  7: 


-32  -  a/1024  -  4(10t r  -  75)  (30  -  2vr  +  16a) 

7  =  - 2(30  -  2.  +  6a) -  (3'58) 


Equation  (3.58)  provides  the  result  for  the  exponential  decay  coefficient  7  for 
an  arbitrary  value  of  the  inclination  angle  a.  These  values  of  7  correspond 
to  the  variational  minimum  of  the  energy  as  a  function  of  a.  The  obtained 
dependence,  being  one  of  our  main  results,  is  presented  in  Fig.3.8. 

Having  obtained  all  the  relevant  numerical  results,  we  are  at  position  to 
discuss  its  physical  properties. 

First,  we  can  conclude  that  a  finite  inclination  angle  a  leads  to  an  increase 
of  both  inverse  normalization  constant  M  due  to  an  extra  finite  region  of 
non-zero  wave  function  (see  Fig.3.1)  and  the  bulk  integral  term  since  another 
piece  of  a  finite  curvature  is  added.  These  regions  were  previously  denoted 
as  i?3  (upper  and  lower).  Apart  from  that,  here  we  are  dealing  with  a  new 
situation,  where  the  variational  minimum  of  the  exponential  decay  coefficient 
7  varies  with  changing  a.  Correspondingly,  in  order  to  apply  the  variational 
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procedure,  we  re-derive  the  value  of  7  for  each  individual  value  of  angle  angle 
a. 

The  variational  results  for  bound  state  energy  are  shown  in  Fig  3.8  at  the 
two  lower  panels.  We  define  the  corresponding  bound  state  energy  as  the 
difference  between  the  obtained  ground  state  energy  and  the  minimal  energy 
of  a  propagating  state.  According  to  our  previous  discussions,  this  critical 
energy  is  equal  to  etreshoid  =  1.234A2/(ma2)  (see  Eq.(2.15)).  We  see  that 
the  ground  state  energy  increases  with  the  inclination  a.  This  increase  is 
finite,  limited  and  does  not  reveal  any  singularities.  Analyzing  our  results,  we 
conclude  that  the  binding  energy  decreases  to  approximately  2/3  of  it  initial 
value.  It  does  not  completely  disappear,  making  the  bound  state  possible  for 
all  values  of  the  inclination  angle.  The  obtained  property  makes  the  S-type 
skew  symmetric  system  distinguished  from  the  standard  X-type. 

3.4  Concluding  remarks 

Finalizing  this  thesis,  I  would  like  to  emphasize  the  uniqueness  and  impor¬ 
tance  of  the  studied  system.  Indeed,  there  are  very  few  structures,  demon¬ 
strating  quantum  bound  states  in  classically  unbound  systems.  In  this  re¬ 
spect  the  problem  definitely  presents  substantial  interest  from  the  theoret¬ 
ical  point  of  view.  Apart  from  that,  simple  two-channel  geometry  enables 
relatively  accessible  fabrication  of  tunable  sensors  based  on  the  electronic, 
transport  and  optical  properties  of  the  S-type  skew  symmetric  cross-wired 
system. 

In  this  thesis,  we  have  analyzed  the  angular  dependence  of  the  lowest 
energy  bound  state  for  an  electron  trapped  at  the  intersection  of  two  identical 
narrow  channels  (quantum  wires)  crossed  at  an  arbitrary  angle  in  such  a  way 
that  the  channel  intersection  area  forms  the  square.  When  the  channels  are 
perpendicular,  this  is  an  unbound  system,  classically,  that  is  known  to  possess 
a  purely  quantum  single  electron  bound  state  [1].  Our  task  was  to  investigate 
the  bound  state  energy  as  a  function  of  the  wire  intersection  angle,  in  order 
to  find  out  whether  the  electron  trapped  state  is  still  there  for  a  non-zero 
inclination  angle,  and  then  to  calculate  the  functional  behavior  of  the  bound 
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state  energy  as  the  channel  intersection  angle  varies. 

The  very  first  main  step  we  have  done  was  to  find  a  right  trial  function 
that  solves  the  problem  by  means  of  the  variational  approach.  Specifically,  an 
analytical  expression  for  the  trial  function  was  found  and  tested.  The  func¬ 
tion  is  supposed  to  be  smooth  in  the  bulk,  and  it  must  drop  down  to  zero  at 
the  boundaries.  In  addition,  the  so-called  roof  lines,  the  lines  of  discontinu¬ 
ous  derivatives,  have  been  found  and  characterized,  and  their  contributions 
have  been  calculated  using  a  special  set  of  recipes. 

Using  the  variational  principle  and  the  trial  function  established,  the 
ground  state  energy  of  the  problem  has  been  calculated.  This  includes: 

•  finding  the  proper  normalization  constant  of  the  accepted  trial  wave 
function, 

•  calculating  the  contributions  to  what  of  the  bulk  integrals  and  those  of 
the  roof  lines.  We  used  the  exponential  attenuation  coefficient  7  as  the 
variational  parameter,  resulting  in  the  lowest  value  of  the  ground  state 
energy  for  any  inclination  angle  a. 

The  bound  state  is  present  for  any  inclination  angle  a.  The  binding  energy 
decreases  with  the  increasing  inclination  angle  a.  However,  the  bound  state 
persists  for  all  accessible  values  of  a,  which  is  known  not  to  be  the  case  for 
other  types  of  the  wire  intersection  geometry.  This  is  the  central  result  of 
our  work. 

Our  results  have  been  compared  with  those  obtained  by  using  other  (sim¬ 
pler)  trial  functions  available  in  quantum  mechanics  texts.  A  qualitative 
agreement  is  demonstrated.  The  results  of  this  thesis  were  presented  with 
success  at  the  annual  American  Physical  Society  meeting  in  Denver,  CO,  in 
March  2014  [49]. 

3.5  Research  outlook 

This  is  mostly  a  theoretical  work,  although  it  maintains  certain  approaches 
and  patterns  of  numerical  experiment,  including  visualizations  and  simula¬ 
tions  using  the  Mathematica  package. 
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The  trial  function  found  was  then  tested  by  using  it  in  numerical  sim¬ 
ulations.  It  was  demonstrated  to  result  in  the  lowest  possible  ground  state 
energy  expectation  value. 

I  would  like  to  summarize  briefly  and  to  discuss  the  importance  of  the 
present  thesis  work,  and  also  to  point  out  its  applications  and  possible  future 
research  directions.  They  are: 

•  This  thesis  work  describes  an  important  example  of  a  quantum  system, 
which  is  unbound  in  a  classical  sense,  but  yet  it  does  possess  a  purely 
quantum  bound  state.  The  results  are  unique  in  that  they  give  the 
angular  dependence  of  the  binding  energy  of  the  electron  trapped  at 
the  skew-symmetric  quantum  channel  intersection.  This  is  something 
that  was  not  reported  before  in  the  literature. 

•  The  data  and  simulations,  presented  in  this  work,  will  be  useful  for 
the  correct  interpretation  of  electron  transport  peculiarities  in  realistic 
quantum  systems  such  as  semiconductor  nanowire  films  and  carbon 
nanotube  bundles. 

•  The  results  obtained  provide  new  useful  information  of  relevance  to  the 
properties  of  semiconductor  quantum  dot-like  systems  and  other  de¬ 
vices  of  confined  geometries.  Understanding  the  properties  of  quantum 
confined  systems  is  important  for  the  development  of  a  variety  of  ap¬ 
plications  in  modern  nanotechnology.  For  example,  a  bound  quantum 
state  in  an  S-type  skew  symmetric  cross-wire  system  could  potentially 
be  used  to  control  electron  transport  on  the  nanoscale,  or  as  an  optical 
sensor,  or  as  a  switch. 

•  Also,  this  thesis  work  presents  an  important  educational  example  that 
is  worth  of  including  into  laboratory  assignments  for  advanced  physics 
courses.  Just  to  mention,  a  simpler  version  of  the  trial  function  for  the 
same  problem  (after  Krishna  Rajagopal)  is  included  in  the  second  edi¬ 
tion  of  the  Introductory  Quantum  Mechanics  textbook  by  D.J. Griffith 
(Pearson,  2005).  My  trial  function  that  has  been  developed  in  this  the¬ 
sis,  has  an  advantage  over  Rajagopal’s  since  it  is  specifically  designed 
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to  minimize  and  simplify  the  roof-line  system,  and  so  it  may  serve  as  a 
pedagogical  example  for  the  critical  development  of  textbook  material 
for  students. 


Appendix  A.  General  approach 
to  roof  lines 


The  mechanism  of  roof-line  integration  is  the  following:  we  find  all  the  exist¬ 
ing  lines,  where  the  wave  function  exhibits  C'^-type  behavior,  which  means 
that  the  function  is  obviously  continuous,  but  the  first  derivative  has  a  finite 
discontinuity  (Heaviside  step).  Consequently,  the  second  derivative  across 
the  roof  line  represents  a  5-function,  thus  must  be  taken  into  account  in  the 
integration. 

We  consider  a  generic  function  in  the  form  of 

F(x)  =  f(x)e~^x)  (4.1) 


and  take  its  second  derivative,  yeilding 


F"(x)  = 


fix)  -  27  f'(x)g'(x)  +  f{x)  (V  ig\x))2  -  7 g'\x ) 


-7 g(x) 


(4.2) 


The  most  general  idea  is  to  find  and  use  only  the  terms,  which  could  poten¬ 
tially  give  a  5  function  and  will  contribute  to  infinitesimally  small  x-integral 

a+e 

/• 

a—e 

Let  us  now  consider  Rajagopal  wave  function  in  the  vicinity  of  the  line 
where  the  exponential  decay  begins.  This  is  the  only  relevant  roof  line  in  our 
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®  a=o 


Figure  4.1:  (Color  online).  Roof  line  schematics. 


case.  For  we  have: 


f(x)  —  1  —  -  9(x  —  a)  +  -9(a  —  x) 
a  L  a 

X 

g(x)  =  9(a  —  x)  H — 9(x  —  a) 
a 


(4.3) 


with  the  second  derivatives: 


f(xY  =  _  a)  > 


(4.4) 


g(x)"  =  -5(x  -  a)  , 
a 

We  collect  only  terms  containing  second  derivatives  of  f{x)  and  g(x): 


d^^x^y) 
dx 2 


=  (f"{x)  ~  7 f(x)g")  e  7  =  — <J(x  -  a) 

CL 


7  (  1  —  —  “ 

a/  a. 


The  result  is  equivalent  to  that  of  the  ”  canonical”  solutions,  provided  in 
the  textbooks.  Obviously,  /  and  g  are  functions  of  both  spatial  variables 
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{x,y},  however  at  this  derivation  we  were  particularly  interested  in  the 
x— dependence. 


- not  a  roof  line 

Figure  4.2:  (Color  online).  Roof  line  schematics. 


For  our  wave  function  the  only  relevant  term  is: 

-7f(x)g"(x)e-^  (4.5) 


g"(x )  =  - 8(x  —  a) 
a 

f(x  a,y)  =Af 

d2'Tb(x  -»•  a,y)  _ 
dx2 


5{x  —  a)e  1 


(4.6) 

(4.7) 

(4.8) 


which  confirms  our  results. 
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